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Abstract We test the feasibility of 3D coronal-loop tracing in stereoscopic EUV 
image pairs, with the ultimate goal of enabling efficient 3D reconstruction of the 
coronal magnetic field that drives flares and coronal mass ejections (CMEs). We 
developed an automated code designed to perform triangulation of coronal loops 
in pairs (or triplets) of EUV images recorded from different perspectives. The 
automated (or blind) stereoscopy code includes three major tasks: (i) automated 
pattern recognition of coronal loops in EUV images, (ii) automated pairing of 
corresponding loop patterns from two different aspect angles, and (iii) stereo¬ 
scopic triangulation of 3D loop coordinates. We perform tests with simulated 
stereoscopic EUV images and quantify the accuracy of all three procedures. In 
addition we test the performance of the blind stereoscopy code as a function 
of the spacecraft-separation angle and as a function of the spatial resolution. 
We also test the sensitivity to magnetic non-potentiality. The automated code 
developed here can be used for analysis of existing Solar TErrestrial RElation- 
ship Observatory (STEREO) data, but primarily serves for a design study of a 
future mission with dedicated diagnostics of non-potential magnetic fields. For 
a pixel size of 0.6"(corresponding to the Solar Dynamics Observatory (SDO) 
Atmospheric Imaging Assembly (AIA) spatial resolution of 1.4"), we find an 
optimum spacecraft-separation angle of a s « 5°. 
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1. Introduction 

One critically needed tool for forecasting severe geomagnetic storms well ahead of 
time is a reliable method to map the magnetic field erupting into the heliosphere 
during coronal mass ejections (CMEs) so that its evolution can be modeled 
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well before their field impacts Earth’s ( e.g. Schrijver et al., 2015). A major 
impediment at present is that we cannot reliably describe the magnetic field in 
nascent CMEs as their erupting structure enters the heliosphere. 

Even with recent advances in the ready availability of vector-magnetic data 
on active regions, the modeling of the configuration of the coronal field above 
these regions remains a challenge. Model results based on surface observations 
alone are generally ambiguous and not well representative of the observed coro¬ 
nal configuration {e.g. DeRosa et al., 2009, and references therein). The use of 
coronal-loop trajectories to guide field models towards a solution compatible with 
surface-field measurements shows promise (Malanushenko et al., 2012, 2014), but 
the fact that only the 2D trajectories projected against the plane of the sky are 
available presents a stumbling block that needs to be overcome. Malanushenko 
et al., (2012) approximate the third coordinate, along the line-of-sight (LOS), 
by pairing up each observed coronal loop with the best-fitting field line of a 
linear force-free held (using separate Held lines for each loop), and they then 
iterate towards a non-linear force-free Held while continually nudging the overall 
solution back to the set of 3D loop trajectories first determined. Whereas this 
results in a model field that by design matches the observed loops quite well, it 
is unlikely that the input 3D trajectories are entirely correct {e.g. see differences 
between 2D and 3D reconstructions in Aschwanden 2013b). Until we have a way 
to measure the 3D loop trajectories we cannot truly validate the method, but 
once the 3D trajectories are known, they can of course be used from the outset 
to guide the model field towards a solution compatible with the observed 3D 
configuration. 

Here, we study a concept with two or three spacecraft that provide stereo¬ 
scopic views of EUV images of coronal loops, which when combined with pho- 
tospheric line-of-sight magnetograms, provide information suitable for 3D re¬ 
construction of the coronal magnetic field. One possible orbital spacecraft con¬ 
figuration is formed by the Sun-Earth Lagrangian points LI and L4 (or L5), 
similar as it was obtained when the Solar TErrestrial RElationship Observatory 
(STEREO) A(head) and B(ehind) spacecraft moved near the L4 and L5 points 
in 2008, while the Solar and Heliospheric Observatory (SOHO) was positioned 
in LI. However, that was a very temporary configuration, with instrumentation 
at moderate resolution. Ideally, a new generation of instruments should have 
a spatial resolution that is comparable to that of Solar Dynamics Observagory 
(SDO)’s Atmospheric Imaging Assembly (AIA), and a field of view at least 
as large as supported by its 4 k x 4 k imaging cameras with 0.6"pixels. Un¬ 
like the STEREO-SOHO configuration, a scientifically more promising mission 
should provide a long-lived multi-spacecraft configuration, providing images with 
essentially identical passband and telescope characteristics. 

We envision that a pair (or perhaps a triplet) of spacecraft equipped with 
the necessary instruments could support the autonomous calculations of the 
coronal magnetic field, its non-potentiality, and the free energies in each active 
region either in near-real time or with up to at most a day delay, so that 
this compound observatory can be used to understand active-region instabilities 
and heliospheric model input, and as an early-warning system for severe space- 
weather storms. To succeed, the data-processing and modeling capabilities would 
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require (i) automated pattern recognition of coronal loops in EUV images, (ii) 
automated stereoscopic pairing of coronal loops, (iii) stereoscopic triangulation 
of coronal loops, and (iv) nonlinear force-free field forward-fitting that yields 
the non-potential magnetic held, its free energy, and (v) - after eruptions - 
quantitative information on the ejecta into the heliosphere. We demonstrate 
the feasibility of the first four of these automated tasks in this study, and we 
constrain the optimum configuration for the angular spacecraft separation. 

Recent reviews on solar stereoscopy and tomography have been presented 
by Aschwanden (2011), and recent reviews on the coronal magnetic held are 
given by Wiegelmann and Sakurai (2012) and Wiegelmann, Thalmann, and 
Solanki (2014). Early attempts of solar stereoscopy using information from a 
single spacecraft (using XUV images from Skylab) used the solar rotation to 
measure stereoscopic parallaxes (Berton and Sakurai 1985), which requires (un¬ 
realistic) static coronal loops on time scales of at least one day, but hydrodynamic 
heating and cooling processes of loops occur on time scales of ss 10 3 seconds in 
active regions ( e.g. Warren and Winebarger 2007), even as the held’s photo- 
spheric boundary is evolving underneath. A dynamic solar-rotation stereoscopy 
method was developed later by Aschwanden et al, (1999, 2000), which relieves 
the requirement of static loops in lieu of a quasi-static magnetic held. This 
assumption is somewhat more realistic, but breaks down after about one day, 
since photospheric magnetic helds involved in major haring and eruptions were 
observed to have characteristic growth and decay timescales of approximately 
one to two days (Pevtsov, Canheld, and Metcalf 1994; Schrijver et al., 2005; 
Welsch, Christe, and McTiernan 2011). Therefore, the only solution for accurate 
stereoscopy requires simultaneous measurements with multiple spacecraft. 

The first stereoscopic reconstruction of coronal loops using two simultane¬ 
ous spacecraft observations was conducted with STEREO/A and B (Feng et 
al., 2007; Aschwanden et al., 2008), but the tracing of coronal loops was carried 
out manually, which is subject to human judgement and does not enable efficient 
processing in real-time, nor in rapid time intervals nor with large statistics. A 
fully automated pattern-recognition algorithm that extracts the 2D geometry of 
coronal loops and performs magnetic modeling has been employed in a recent 
study (Aschwanden, Xu, and Jing 2014), applied to 172 hare events in numerous 
active regions. This algorithm also performed nonlinear force-free held modeling 
and determined the evolution of the free energy during hare events, using high- 
resolution images of SDO/AIA, but this algorithm uses information on the 
2D geometry from a single spacecraft only, and thus is expected to retrieve 
less accurate information on the magnetic held than would be possible from 
stereoscopically determined 3D geometries of coronal loops (Aschwanden 2013b). 
Comparisons of NLFFF reconstructions using single-spacecraft 2D versus dual¬ 
spacecraft 3D geometries yielded consistent results for a simple forward-htted 
quasi-NLFFF model in terms of vertical currents (Aschwanden 2013), but the 
accuracy for more general NLFFF solutions is not known. On the other hand, 
accurate NLFFF solutions htted to coronal loop geometries have been accom¬ 
plished with a Quasi-Grad-Rubin method (Malanushenko et al., 2012, 2014), 
but the htting constraints were based on manual tracing of loops and the com¬ 
putation time with the present code prevents efficient real-time calculations. 
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Even observations from STEREO/A and B during the most optimum con¬ 
ditions at small spacecraft-separation angles (during 2007) were not able to 
provide accurate magnetic-loop geometries, because the spatial resolution of 
the STEREO Extreme UltraViolet Imager (EUVI) images is too poor, being 
three times poorer than SDO/AIA images. Given all of these instrumental and 
computational restrictions, an ideal multi-spacecraft configuration suitable for 
most accurate magnetic-field modeling and optimum signal-to-noise ratio has 
still to be established with a new design for a future mission. 

In this article we test the principle of dual and triple stereoscopy to establish 
the 3D coronal loop configuration from simultaneous EUV image combinations. 
In the process, we develop a suite of numerical codes that is capable of perform¬ 
ing stereoscopy in an automated way, which includes simulations of synthetic 
stereoscopic image pairs (section 2), automated detection of coronal loops in 
high-resolution EUV images (section 3), automated stereoscopic pairing of loops 
(section 4), and automated 3D triangulation of loops (section 5). We investigate 
the accuracy of stereoscopy as a function of the number of spacecraft (section 
6), as a function of the spacecraft-separation angle (section 7), as a function 
of the spatial resolution (section 8), spacecraft position (section 9), and its 
sensitivity to the non-potentiality of the magnetic field (section 9). Discussions 
and Conclusions are provided in section 10. 


2. Simulation of Stereoscopic Images 

For our simulations of EUV images suitable for testing the principle of multi¬ 
spacecraft stereoscopy we choose data from active region NOAA 11158, as ob¬ 
served on 15 February 2011, around the time of a GOES X2.2-class flare event 
that occurred at 01:56 UT. This active region produced the first X-class flare 
event in the era of the SDO (Pesnell, Thompson, and Chamberlin 2011), and 
this is one of the best studied regions. This active region was also chosen in 
previous magnetic modeling with nonlinear force-free field (NLFFF) methods 
(Malanushcnko et al, 2014), using Helioseismic and Magnetic Imager (HMI) 
(Scherrer et al, 2012) and AIA (Lemen et al, 2012) data. 

In our simulation of image sets from different perspectives, we start with 
a line-of-sight magnetogram of HMI only, acquired at 15 February 2011, 01:40 
UT. Observed from Earth perspective, the center of AR 11158 has a heliographic 
position of S21W12, which is centered at [xo,yo] ~ [0.10, —0.35]R e from disk 
center. The HMI image has a pixel size of 0.50422" and the solar radius is 
1927.2 pixels. We extract a subimage with a field-of-view x = [0.0,0.5]Rq and 
y = [— 0.5, 0.0]Rq, which corresponds to a size of 965 x 965 HMI pixels. 

We decompose the magnetogram into n mag = 100 sources and calculate the 
potential field that results from the combined field of the 100 subphotospheric 
unipolar magnetic charges, according to the method described in Aschwanden 
and Sandman (2010) and in Appendix A of Aschwanden et al., (2012b). The line- 
of-sight magnetic field component has a range of —700 G < B z (x, y) < 1017 G 
within the chosen field-of-view. We cover the 965 x 965 pixel subimage with a grid 
of 100 x 100, and define each grid point that has a LOS field B z (x , y) > 200 G as 
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Figure 1. Simulated pair of stereoscopic images E15 and E00, seen from an angle 15° East 
from Earth (right bottom panel) and from Earth (left bottom panel). A heliographic coordinate 
grid is indicated with increments of 5°. The simulated EUV images are composed from 261 
magnetic-field lines (top panels), generated from a potential-field extrapolation of a magne¬ 
togram of active region NOAA 11158 observed with HMI/SDO on 15 February 2011, 01:14 UT. 
Only field lines with magnetic field strengths of \B Z \ > 200 G at the footpoints are displayed. 


a footpoint of a coronal loop, from which we extrapolate the potential-field line 
until it hits one of the six boundary sides of the computation box, using a height 
of /imax = 0.15 Rq (or 104 Mm). From the 10,000 grid points, the magnetic field 
exceeds the minimum limit of 200 G at 261 locations, which yields 261 potential- 
field lines. We rotate the computation box to different viewing angles, using the 
same coordinate transformation as solar rotation produces, for instance rotating 
by +15° to the West, in order to mimic a viewing position of STEREO-B at 
position E15 eastward on the Earth (Figure 1, top panels). 
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The barometric density is n e (h) oc exp(— h/X), with the scale height A = 50 
Mm x T e [MK], corresponding to a temperature of T e = 2 MK, which is typical 
for structures that are seen in the AIA 193 A and 211 A images. The intensity 
of the image scales with the emission measure, i.e. F{x,y) oc f EM dz oc n^Az, 
with A z the LOS-integrated column depth. For sake of simplicity, we do not 
intend to simulate AIA images at particular wavelengths, because the results of 
stereoscopic simulations depend primarily on the geometry and signal-to-noise 
ratio of the detected loops, which applies to any temperature or EUV wavelength. 

In order to create an adequately realistic EUV image, we convolve each point 
of a field line with a Gaussian kernel that mimics typical loop aspect ratios (of 
the loop width to the length) and gravitational stratification. For the half width 
of a loop we choose the scaling of w(s) = wq y/n^, where n s is the length of the 
loop in pixels, and wq = 1 Mm is the minimum loop width. In order to mimic the 
point spread function (PSF) of the instrument, we simulated loop widths [w(s)] 
with Gaussian kernels that are always larger than the PSF, i.e. w(s) > w ps { ~ 
two pixels. 

Simulated EUV flux maps of the optically thin plasma are rendered in Figure 
1 (bottom panels), similar to the method of Gary (1997). In the later sections 
of this article we also simulate similar EUV maps with different spacecraft- 
separation angles (section 6 and 7), with different spatial resolutions (section 8), 
or with different magnetic field models (section 10). 


3. Automated Pattern Recognition 

Stereoscopy of coronal structures has been pioneered only by visual tracing so far, 
for instance using Skylab images (Berton and Sakurai 1985) or STEREO/EUVI 
image pairs (Feng et al., 2007; Aschwanden et al. 2008). Manual tracing of coro¬ 
nal loops, however, is very time-consuming and depends on human judgement, 
and thus is very inefficient and subjective, preventing any frequent sampling or 
real-time operation. The availability of an automated pattern-recognition code is 
therefore a valuable element to accelerate progress in magnetic-field modeling of 
the solar corona. In our context here, we aim for a “blind stereoscopy method”, 
rather than “stereoscopy aided by visual guidance”. 

Five experimental numerical codes for automated tracing of coronal loops 
were compared in an initial study (Aschwanden et al., 1998). One of them, 
the so-called oriented coronal curved loop tracing (OCCULT-1) code was fur¬ 
ther developed and approached visual perception (Aschwanden, 2010). A new 
advanced code (OCCULT-2) was further optimized for curvi-linear tracing ap¬ 
plied to Transition Region and Coronal Explorer (TRACE) data (Handy et 
al, 1999), SDO/AIA data, Swedish Solar Telescope (SST) data, and microscopic 
biophysics images (Aschwanden, DePontieu, and Katrukha 2013). Here, we use 
this OCCULT-2 code. The automated pattern recognition algorithm detects 
iteratively curvi-linear patterns with large curvature radii, starting at a position 
with the highest flux, propagating along the local ridge guided by the local 
curvature radius, and it erases the signal of a traced loop segment from the 
image before it starts with the next loop segment. The automated loop tracings 
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Figure 2. Automated loop tracing with the OCCULT-2 code (red curves), superimposed on 
the simulated 261 magnetic-field lines (black curves in top panels) and the simulated EUV 
images (bottom panels). A total of 80 (or 84) loop segments are detected at E00 (or E15) 
above a flux threshold of 1% and with a minimum length of / m ; n > 20 pixels. 


are carried out here independently in each image of a stereoscopic pair (such as 
E00, E15 shown in Figure 2). 

For an example, we show the automated tracing of a pair of simulated stereo¬ 
scopic images in Figure 2, which was simulated using 261 magnetic field lines. 
The OCCULT-2 code detects from the EUV images a total of 80 loop segments 
with a length of l m in > 20 pixels in image E00, and 84 loop segments in image 
E15. The parameters can be adjusted depending on the type of data or desired 
pattern. What is particular about the simulated stereoscopic images here (Figure 
1) is freedom from noise, in contrast to observed EUV images. Noise-free images 
allow for more sensitive detection of faint structures and are less prone to mis- 
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guided detections in faint structures that are comparable with the ambient noise 
level. We simulate noise-free images here in order to study the performance of 
automated stereoscopy under ideal conditions, but will add data noise later to 
study the stereoscopic behavior under more realistic conditions. 

For our application here we chose the following parameter settings for the code 
OCCULT-2: a pixel size of 1.5" (or 1.0 Mm), a highpass filter of nsm\ = 1 pixel, 
a minimum curvature radius of r m j n = Z m ; n pixels, a minimum loop segment 
length of Z m i n = 20 pixels, an image base level of thr\ = 0, a threshold level of 
thr 2 = 0.01 (in units of the highpass-filtered flux maximum), no gap (n gap = 0) 
along a traced structure, and a maximum loop number of n max = 200 loop struc¬ 
tures per image. These settings yield a near complete detection of unconfused 
loop segments, down to the faintest structures seen visually (Figure 2). Since 
structures are generally seen down to the spatial resolution of the instrument, 
the same settings in units of pixels are recommended also for an instrument 
with a different spatial resolution (although it corresponds to a different abso¬ 
lute scale of the pixel size.) The most challenging part is the crowded central 
core of the active region, where multiple loops overlap and cross each other. A 
correct disentangling of loop structures in such nested areas can probably only be 
achieved by forward-fitting of multiple loop geometries, rather than by iterative 
loop tracing. 


4. Automated Stereoscopic Pairing 

The second major task of the autonomous stereoscopy procedure is the pairing 
of corresponding loops, namely the correct association of loop segment i in 
image A with the stereoscopic counterpart of loop segment j in the stereoscopic 
image B. This problem of “stereoscopic correspondence” or “stereoscopic pairing 
ambiguity” has never been systematically investigated, and thus we explore it 
here to some degree to enable automated stereoscopy. 

Stereoscopy is generally accomplished by transforming a stereoscopic image 
pair into an epipolar coordinate system (Inhester 2006), which generally requires 
a rotation and rescaling of each image (if the images are taken with different 
image scales or at different distances from the Sun). The epipolar plane is 
defined by three points: the Sun center, and the positions of two stereoscopic 
vantage points, which are given by the spacecraft locations A and B for the solar 
STEREO mission. In such an epipolar coordinate system, an image A is taken in 
the [x, y]-plane, the epipolar rotation axis is in the y-direction, which warrants 
that the stereoscopic parallax causes a rotational shift in the ^-direction only, 
while the y-coordinate remains unchanged. Consequently, any structure with 
coordinate [*(s),y(s)] in image A, where s is a loop-length coordinate, has the 
coordinates [x(s) + Ax(a s ) + A x(h),y(s)] in image B, where Ax(a s ) is the 
rotational shift of the coordinate system due to the spacecraft-separation angle 
a s , and Ax[/i] is the parallax that depends on the altitude h above the solar 
surface. The range of possible altitudes, 0 < h < /i max , defines the solution space 
in image B, where a corresponding loop can be located. In the ideal case, a loop 
is detected over its entire length in both images, and appears isolated within 
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the search area. The search area for such a loop with coordinates [xi(s), yi(s)] 
in image A, is bound by [^(s), 2 / 2 (s)] with yi(s) = j/ 2 (s) in the y-direction and 
£i(s) < |iC 2 (s) — A:c(a s )| < Aa;(ft max ) in the x-direction. If there is only one 
loop segment in this search area in image B, the correspondence is unique and 
a stereoscopic triangulation can directly be calculated. In reality, however, there 
are often multiple loops in the search area and we have to develop a strategy to 
find the most likely stereoscopically correct correspondence. 

We illustrate the “stereoscopic correspondence problem” in Figure 3, where 
we show stereoscopy between a spacecraft W15 and Earth view (E00) (Figure 
3 left panels), as well as stereoscopy between a spacecraft E15 and Earth view 
(E00) (Figure 3 right panels). A total of n = 80 loops were detected in image 
(E00), while 84 loop segments were detected in image (E15), whereof nc = 
51 segments have an overlapping y-range in both images. The tracing of the 
largest loop in image E00, which has the number #16 when sorted by length, 
is outlined (black/white dashed linestyle) in Figure 3 in all panels. Rotating a 
loop structure from image E15 onto the view of E00, using two fixed distances 
from Sun center (i.e. with a minimum altitude ft = 0 and a maximum altitude 
of ft max = 0.15R q ), we find a solution space in E00 for each structure detected 
in E15. In this case we find two loop segments that overlap with the traced 
loop #16 in image E00, so there is a two-fold ambiguity which loop should be 
stereoscopically triangulated. The boundaries of the solution space of loop #35 
from image E15 is indicated in the image E00 with a red-blue zone, where blue 
corresponds to a minimum altitude of ft = 0 and red to a maximum altitude 
of ft = ftmax = 0.15 solar radii. Which is the correct correspondence? Since the 
probability of a true correspondence increases with the length of the coincident 
segment, we use the criterion of maximum length, which indeed corresponds to 
the correct solution (indicated with an orange line in the top panels), known 
from the simulated magnetic-field lines. 

How large is the ambiguity of pairing stereoscopic loop segments? We count 
the number of loop segments in E15 to each loop segment of E00 that intersects 
with the stereoscopic solution space, bound by an overlapping y-range and alti¬ 
tude range ft = [0, ft max ]. We find that most loops have an ambiguous stereoscopic 
correspondence, within a range of 1-10 possible correspondences, or a statistical 
mean of n am b = 3.2 ± 2.3. In the example shown in Figure 3 (right panel), 
there are two ambiguous loop segments in image E15 that could potentially 
correspond to loop #16 in image E00 within the altitude range used. The degree 
of ambiguity generally depends on the specified altitude range, which is chosen 
to be ft m ax = 0.15 solar radii here, and is expected to linearly increase with larger 
altitude ranges. One strategy to reduce the number of ambiguities is to eliminate 
those that have already been used previously in the iterative stereoscopic pairing. 
A further strategy to avoid false stereoscopic pairings is to start with those that 
have the longest loop segments, where the least ambiguity occurs. Proceeding 
to smaller and smaller loop segments, the number of ambiguities then decreases 
systematically. 

Based on these considerations, we implement the following steps in the (blind) 
stereoscopy code as a strategy to optimize the stereoscopic pairing procedure. 

(i) All riA detected loops in A are sorted by their length. 
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Figure 3. The largest loop traced in image E00 (second panel from bottom) is #16 
(black/white dashed linestyle), which corresponds to loop #9 in image W15 (bottom left 
panel), or to loop #35 in image E15 (bottom right panel). The projection of loops rotated at 
a photospheric level h = 0 is indicated with blue curves, and rotated at a (maximum) coronal 
height of h = O.15R.0 with red curves. The side views (top panels) indicate the stereoscopically 
triangulated altitudes of candidate loops in the solution space. The correct solution known from 
the simulated magnetic field is indicated with a field line drawn in orange color. 
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(ii) All ns detected loops in B are sorted by their length. 

(iii) For each sorted loop i\ = 1, n\ we determine which of the sorted loops 
*b = 1, ...,n b overlap with the solution space of loop within an altitude 
range of h = [0,h max ]. 

(iv) The longest segment in B that overlaps with the altitude range of loop i\ 
is selected as the stereoscopic counterpart *B,sei- 

(v) A loop *B,sei in B that has already previously been paired with a loop *a, 
is excluded for pairing in the next iterative pairing step. 

An alternative approach to solve the loop correspondence problem is the 
so-called “magnetic stereoscopy method” (Wiegelmann et al., 2006; Feng et 
al., 2007), where an extrapolated linear force-free (LFFF) magnetic held is used 
to identify corresponding loops. A possible advantage of this method is that 
it reduces the solution space of corresponding loop locations more efficiently 
than our empirical stereoscopic pairing method described above, especially for 
large spacecraft-separation angles. However, a disadvantage of this method is 
that the initially chosen magnetic held (LFFF) model introduces a bias that 
favors solutions close to the initial (LFFF) model and may even prevent the 
convergence towards a nonlinear force-free held (NLFFF) solution. However, the 
best method may be an iterative approach, where stereoscopic loop pairing and 
magnetic held modeling is performed in alternating steps, starting from an initial 
potential held model, a and ideally ending at a best-htting hnal non-potential 
held model. 


5. 3D-TYiangulation of Loops 

The third step in the blind-stereoscopy procedure consists of the triangulation 
of loop points. This is the easiest part of the stereoscopy procedure, because it is 
a uniquely dehned mathematical geometry problem, after we have identihed 
the correct corresponding loop counterparts in both images A and B in an 
epipolar coordinate system. Specihc triangulation formula are given in a number 
of previous studies ( e.g. Berton and Sakurai 1985; Inhester 2006; Aschwanden 
et al, 2008). Geometric parameters of the heliographic coordinate system of 
an image are generally specihed in the FITS headers of the image data hies 
(Thompson and Wei 2010). 

Here, we derive the analytical relationships in their simplest form for a pair of 
two images A and B that have been already rotated and scaled into an epipolar 
coordinate system, which contains only four variables for every loop point s: 
(x\ ,y) are the coordinates of a loop point in image A with respect to the Sun 
center, (sb, y) are the coordinates of the corresponding loop point in image B, 
where the y-coordinate is identical in an epipolar coordinate system (y = ua = 
Ub), and the stereoscopic spacecraft-separation angle [a s ], measured from Sun 
center in the epipolar plane. Scaling the distances in units of solar radii, the 
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distance from Sun center is r = 1 + /i, where h is the altitude above the solar 
surface. The distance of point (xa ,y) from the solar (epipolar) axis is then 

p = cos b (1 + h) , (1) 


where ( l,b ) are the heliographic longitude and latitude in a Stony hurst grid 
(with l = 0 and b = 0 at solar disk center). The cartesian coordinates xa,2/a, 
Xb are then related to the heliographic coordinates 1 aJb,& by the following 
relationships, 


x A = p sin (1 a) , 

(2) 

x B = p sin (l B ) , 

( 3 ) 

Ib = Ia + 0! s , 

( 4 ) 

l = sin (6) (1 + h) . 

( 5 ) 


Now we can substitute and eliminate the observables (1a, 1b, 2/ ) and obtain the 
relationships for the variables ( p 1 b, h ), 


P = 


x\ 


xb — xa cosa s 
sina s 


1/2 


( 6 ) 


b = arctan 



( 7 ) 


h = 


cos b 


- 1 . 


( 8 ) 


Another simple method is to rotate the coordinate of a location (xb , y ) from 
image B by the spacecraft angle a s into the coordinate system of image A 
using two different altitudes \h\ and I 12 ], so that the correct altitude h can 
be interpolated at the matching position x'g* = xa- We used both methods in 
order to validate our triangulation code. 

A result of stereoscopic triangulation of loops measured in the images E00 and 
E15 is shown in Figure 4, along with their projections into orthogonal planes. 
The 80 loop segments that were automatically traced in image E00 are displayed 
in Figure 4c (thin solid curves), while those segments for which a corresponding 
match in E15 was found are indicated with thick solid curves. The triangulated 
heights are shown as a function of the x-coordinate, r[x) = 1 + h(x) (Figure 
4b), and as a function of the y-coordinate, r(y) = 1 + h(y) (Figure 4d). The 
closest matching 80 magnetic field lines that match the loops traced in image 
E00 are indicated also (Figure 4a and 4e). The display in Figure 4 demonstrates a 
good match between the theoretical magnetic field lines and the stereoscopically 
triangulated loops. 

The numerical accuracy of stereoscopic triangulation depends somewhat on 
the spatial direction of the coronal loop or magnetic field line. In the epipolar co¬ 
ordinate system, the stereoscopic parallax occurs in the x-direction, which yields 
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E00= 80, W15= 85, Stereo= 52, Length > 20 pix 
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Figure 4. Automatically traced loops in image E00 (thin colored curves) and loop segments 
with stereoscopic correspondences in E00 and E15 (thick colored curves) are shown in the x — y 
plane (panel c). The stereoscopically triangulated solutions are shown in three projections, in 
the x — r plane (panel b), and in the r — y plane (panel d). The theoretical solutions known 
from the magnetic field simulations are shown in the panels a and e. 


the most accurate measurement if a loop or field line is oriented in the y-direction, 
i.e. in the North-South direction. If the loop is oriented in the x-direction, there 
is a singularity in the altitude inversion, because the parallax direction coincides 
with the loop direction, and the correspondence of a loop segment in a pair of two 
stereoscopic images is mathematically ill-defined, which prohibits stereoscopic 
triangulation at this location. This singularity, which we may call the “epipolar 
degeneracy”, can affect the accuracy of stereoscopic triangulation for a range 
of angles where the tilt-angle tan(|$(s)|) = |y(s i+ i) — y(sj)|/|x(si + i) — a:(si)| 
along a loop coordinate s has a small value, due to the finite spatial resolution 
and directional tracing errors (Aschwanden et al, 2008, 2012b). In order to 
overcome this epipolar degeneracy problem, we apply stereoscopic triangulation 
only at loop locations [s] where the loop direction is larger than a critical value, 
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i.e., |tan$(s)| > 0.1, and apply a low-order polynomial interpolation for the 
coordinate z(s) in those gaps. 


6. Dual Versus Triple Spacecraft Stereoscopy 

The minimum option for solar stereoscopy is two perspectives, but one may 
opt for a third for a number of reasons. First of all, should one instrument 
fail, one has still a fall-back option with two spacecraft that can essentially 
accomplish the stereoscopy task. Secondly, stereoscopy with two perspectives 
often confronts us with the problem of ambiguous correspondence. Which loop 
structure from perspective A has to be triangulated with what loop structure 
from perspective B? The combination of three perspectives virtually eliminates 
the stereoscopic-ambiguity problem. 

We demonstrate the bootstrapping effect of stereoscopic triangulation as it 
could be achieved with three spacecraft in the example shown in Figure 5. The 
view of the active region from the three spacecraft E00, E15, and W15 is depicted 
in Figure 3, and stereoscopy of loop #16 from the image pair E00 and E15 
(Figure 3 right panels), and the pair E00 and W15 (Figure 3 left panels) yields 
a self-consistent solution that is close to the theoretical values (orange curves in 
Figure 3). In Figure 5 (top left panels) we show the orthogonal projections h(x) 
and h(y) for the same loop #16, found from the spacecraft pair E15+E00 (Figure 
5, blue curves), and from the spacecraft pair W15+E00 (Figure 5, red curves), 
along with the theoretical solution (Figure 5, orange curves in dashed linestyle), 
which all agree within Ah < 0.01 solar radii. In the same representation we show 
the same information for the 12 longest detected loops in Figure 5. From these 
12 stereoscopic triangulations we see a consistent solution of both spacecraft 
pairs with the theoretical model field lines in eight cases, while the results from 
E15+E00 (Figure 5, blue curves) fail for the four loops #18, 19, 31, and 45. 
Nevertheless, the results from W15+E00 are correct in 11 out of 12 cases. For 
this particular example, which may be typical for many other observations, we 
can say that stereoscopic triangulation with three spacecraft is successful in 
« 90% (11 out of 12 cases), while triangulation with two spacecraft can have 
a reduced success rate in the range of 75%-90%, depending on the position 
of the spacecraft. A generalization of our two-spacecraft stereoscopy code to a 
triple-spacecraft configuration could easily be implemented, based on the relative 
overlap range of the y(s)-values of the automatically traced loop segments in each 
of the three spacecraft images, and the highest probability or correct stereoscopic 
correspondence based on the maximum lengths of the paired loop segments 
among the three images from different vantage points. 


7. Spacecraft Separation Angle 

What is the optimum spacecraft-separation angle for stereoscopy? In a previous 
study with STEREO/EUVI data, it was demonstrated that stereoscopic trian¬ 
gulation is in principle possible from small (6°) to large (170°) angles, based on 
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Figure 5. Comparison of solutions for altitudes as a function of the ^-coordinate, h(x ), and 
the y-coordinate, y(h), from the stereoscopic triangulation between the images E00 and E15 
(blue curves), and the images E00 and W15 (red curves). The theoretical solutions based on 
the closest simulated magnetic-field lines are indicated in orange. Inconsistent solutions occur 
for four loops (#18, 19, 31, 45) out of the 12 cases, but are correct for the spacecraft image 
W15 (red curves) in all cases but #19. 
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a small sample of visually traced loops (Aschwanden et al., 2012b). Combining 
the accuracy of altitude triangulation with the stereoscopic correspondence am¬ 
biguity, it was estimated that a spacecraft-separation angle of a s = 22° — 125° is 
most favorable for stereoscopy, using an instrument with the spatial resolution 
2.6" (with pixel size of 1.6") such as STEREO/EUVI. 

Here, we simulate data for a spacecraft-separation angle of a s = 1° to 90° in 
both the eastern and western direction, and perform stereoscopic triangulation 
between a spacecraft at angle a s (W90,...,W01, E01, ...E90) and the spacecraft at 
Earth view (E00). For simplicity we label the positions at angles (W90,...,W01, 
E01, ...E90) with “spacecraft A”, and the position at Earth view with “spacecraft 
E”. 

For each of the spacecraft positions we measure the number of automatically 
traced loops (Figure 6a), which has a maximum value of Nq = 80 (in image E) 
and tends to decrease as the region is seen closer to the solar limb (in the images 
A). If the distribution of loops were homogeneous and isotropic in an active 
region, the view would be rotation-invariant and the number of detected loops 
should be constant as a function of the aspect angle. Therefore, the observed 
slight decrease of detected loops towards the limb indicates a larger horizontal 
than vertical extent of the active region. If we assume a homogeneous density 
of loops in a box with a horizontal length [Ax] and height [h max ], the projected 
length of the box as a function of the rotation angle (or longitude) a s , we expect 
that the number of detected loops is roughly proportional to the projected length 
of the box, for which we then expect a center-to-limb variation of 

„ , v - r Axcos(a s ) + h max sm(a s ) 

Ni oop {a s ) « N 0 - — - , (9) 

where No is the maximum number of detected loops, Ax is the East-West exten¬ 
sion of the active region, and h max is the maximum altitude. We overplot such 
a function in Figure 6a, using Ax = 0.5Rq based on the chosen field-of-view, 
and h max « 0.15Rq, which approximately reproduces the decrease of detected 
loops towards the limb. The AR is located at 12° West, causing an obscuration 
by the limb at —78° and the number of detected loops to go to zero beyond the 
limb. 

Then we measure the triangulation efficiency qt T i in suitable loop segments 
(Figure 6b). This number is defined by the ratio qt T i = -Wri/A^aii of the number 
of loop positions N tl i where a valid stereoscopic triangulation could be executed, 
normalized to the total number 7V a n of all possible loop positions. The require¬ 
ment for a stereoscopic triangulation of a loop position is a valid stereoscopic 
correspondence between two spacecraft A and E, which is an identical y-position 
in both spacecraft A and E, and a valid altitude range of 0 < h < /i max in the 
stereoscopic triangulation. For this number we find a typical value of qtri ~ 0.5, 
mostly caused by incomplete loop tracing or erroneous stereoscopic correspon¬ 
dence identifications. As we can see in Figure 6b, the number of triangulated 
loop positions qtri(ct s ) follows a similar function as the number of detected loops 
(Equation 9), with a drop-off at both the eastern and western limb. Thus the 
efficiency of stereoscopy is warranted in a broad range of |a s — Zq| ~ 60° between 
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Figure 6. (a) Statistics of the number of automatically traced loops (top panel), (b) the 
number of triangulation points (second panel), (c) the mean number of ambiguous loops per 
triangulated loop, given by the mean and standard deviation (third panel), and (d) the median 
accuracy of altitude measurements dh in units of solar radii, as a function of the spacecraft 
heliographic position (or spacecraft-separation angle from Earth). The colored curves corre¬ 
spond to the theoretical model described in the text. The spatial resolution corresponds to a 
pixel size of Ax = 1.4 /; . 
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Figure 7. The number of ambiguous loops in stereoscopic correlation as a function of the 
spacecraft-separation angle [a s ] is shown. If spacecraft A located at a = 0° detects one single 
loop within a spatial width of Ax = 1, a spacecraft B at position a = 90° detects a number 
of jV am h (a) = VN loops that all have the same projection after they have been rotated by an 
angle of —a to the viewpoint of spacecraft A in the ^-direction. For an intermediate position 
B at q, the number of ambiguous loops scales as Af am b(a) = %/~N sin(a). 


a spacecraft position A and E, where Iq ~ 12° is the longitude of the active 
region. 

The number of ambiguous loops may also play a significant role in the evalu¬ 
ation of the optimum spacecraft angle for stereoscopy. We performed automated 
stereoscopy between a near-Earth spacecraft E and a spacecraft A at any position 
from the most western viewpoint at a s = —90° to the most eastern view¬ 
point a s = +90°, in increments of 5° for the whole range, and in increments 
of 1° in the small-angle range of —10° < a s < +10°. While the automated 
code detected N\ oop = 80 loop structures in image E, an average number of 
1 < nioop.A ~ 10 ambiguous loops were detected in image A, where “ambiguous” 
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Figure 8. The error trapezoid of stereoscopic triangulation is shown (grey area), given by 
the two lines of sight of the two observer directions A and B, separated by an angle a s . The 
uncertainties [An] in the n-direction and in the z-direction depend on the pixel width [A p ; x ] 
and half aspect angle a. s /2 (Aschwanden et al., 2012b). 


means the number of candidate loops in image A that have a valid altitude range 
0 < h < h max = 0.15 in the stereoscopic triangulation. The variation of the 
number of ambiguous loops is shown in Figure 6c, which reveals a minimum of 
one single loop at a s ss 0°, while the ambiguity seems to increase up to angles of 
a s ~ 50° — 70°. We can model the number of ambiguities by assuming a uniform 
distribution of North-South oriented loops in a cube, which is most favorable 
for stereoscopy (Figure 7). If we have n x loops in ^-direction and nh loops in h- 
direction, the total number of loops is Ni oop = n x xrih- For a quadratic box with 
a maximum number N\ oop of detected loops we have n x = rih = \JN \ oop . The 
number of ambiguous loops scales then with the projected length (Ax in Figure 
7) as a function of the rotation angle (or spacecraft-separation angle a s ). If a 
loop were detected within a width of Ax = 1 (Figure 7) for a small stereoscopic 
viewing angle, the projected width of all stereoscopically corresponding loops 
is Ax « VN at the limb (Figure 7), and scales according to the sine-function 
inbetween, 

rfam6(rfs) — IT \J'Nioop sin|cy s | , (10) 

which is overplotted on the measurements in Figure 6c. The predicted number 
of ambiguous loops matches the numerically determined number fairly accurate 
over the range of —50° < a s < 80°. The ambiguity factor is almost symmet¬ 
rical for a western and eastern separation angle, but exact symmetry is not 
expected for an active region that has no symmetry in its magnetic field and 
EUV brightness. 

The ultimate parameter that determines the optimum spacecraft-separation 
angle for for dual-spacecraft stereoscopy is the accuracy of altitude measure¬ 
ments, which depends not only on the ambiguity factor (for stereoscopic cor- 
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respondence) but also on the spatial resolution of the images. If we make the 
spacecraft-separation angle smaller and smaller, the horizontal parallax as a 
function of the altitude becomes also smaller and will reach sub-pixel scale, 
where the stereoscopic information becomes unmeasurable, because there is a 
singularity when the spacecraft-separation angles approach zero. Theoretically, 
the error in the line-of-sight extent [Az] of a point source observed with a 
spatial resolution [Ax] can be understood from the “trapezoid relationship” 
(Aschwanden et al., 2012b) shown in Figure 8, 


A h res — 


Ax/2 


sin (|a s |/2) 


( 11 ) 


which exhibits a singularity at spacecraft-separation angle a s = 0 in the number 
of stereoscopic triangulation points (Figure 6b) and in the accuracy [A/i] (Figure 
6d) of stereoscopically triangulated altitudes. For larger spacecraft-separation 
angles, this finite spatial resolution effect becomes negligible, while uncertainties 
due to the ambiguity factor in stereoscopic correspondence dominate. Thus, the 
center-to-limb variation of the accuracy of altitude measurements is expected 
to vary as a sine-function as the ambiguity function does (Equation 10), where 
the maximum (positive or negative) error corresponds to the half height range 

(^max/2), 

A z amb = fhtAA sin (|a s |) . (12) 

In order to obtain an error in altitude [Ah], we have to correct for the cosine-angle 
of the projection between the line-of-sight z and the altitude h, 


A h arn i, — A z arn b cos |a s | , 


(13) 


which yields the combined error (added in quadrature), 


Ah = 


A hi. 


Ah l mb = 


( Ax/2) \ 2 

\sin(KI/2)y 



a s cos a s 


2 


1/2 


(14) 

This theoretical prediction is overplotted on the numerical datapoints of the 
uncertainties for the stereoscopic altitude measurements (Figure 6d, red smooth 
curve), which matches the data closely in the range of —60° < a s < +60°. An 
interesting consequence of this model is that it predicts a highest accuracy at 
a spacecraft separation angle of a s ~ 8.2°. According to the simulations, the 
highest accuracy is Ah ~ O.O2R 0 (or 14 Mm), does not deteriorate more than 
about a factor of two to Ah < O.O4R 0 (or 28 Mm) at larger separation angles. 

This is also the reason why a relatively wide range of stereoscopic angles 
of a s « 22° — 125° was found to be usable for stereoscopy according to an 
earlier study (Aschwanden et al., 2012b), a range that has been adopted for 
another proposed future stereoscopic mission (Strugarek et al., 2015). Note that 
the uncertainty of the ambiguity was estimated differently in the previous study, 
assuming a dependence of A/i am b oc 1/| cos (a s )| based on a theoretical argument 
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Figure 9. Same representation as Figure 6, but for a spatial resolution corresponding to a 
pixel size of Ax = 0.6"): (a) Statistics of the number of automatically traced loops (top panel), 
(b) the number of triangulation points (second panel), (c) the mean number of ambiguous loops 
per triangulated loop, given by the mean and standard deviation (third panel), and (d) the 
median accuracy of altitude measurements dh in units of solar radii, as a function of the 
spacecraft heliographic position (or spacecraft-separation angle from Earth). 
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Table 1. Dependence of optimum spacecraft-separation angle [cKbest] and the uncer¬ 
tainty of stereoscopic altitudes [Ah] on the spatial resolution [Ax] and altitude range 
[hmax] of the solution space. 


Instrument 

Pixel 

Spatial 

Pixel 

Altitude 

Spacecraft 

Altitude 


size 

resolution 

size 

range 

separation 

uncertainty 






angle 



Ax 

PSF 

PSF 

hmax 

CX-best 

Ah 


[arcsec] 

[arcsec] 

[Mm] 

[Mm] 

[deg] 

[Mm] 

AIA 

0.6 

1.4 

0.44 

70 

05 

CO 

o 

5.5 

EUVI 

1.6 

2.6 

1.16 

70 

10.7° 

8.9 

AIA 

0.6 

1.4 

0.44 

140 

4.6° 

7.8 

EUVI 

1.6 

2.6 

1.16 

140 

7.4° 

12.6 


(Equation 7 in Aschwanden et al., 2012b), while we obtain a probably more 
realistic dependence of A/i am b oc sin (a s ) cos (a s ) (Equations 12 and 13) here, 
based on numerical simulations of stereoscopic triangulations. The large range 
of a s = 22° —125° still specifies an angular range where stereoscopy is feasible, if 
we tolerate a factor of two in the uncertainty of altitude measurements (Figure 
6d and 9d), but the optimum angle of 4.6° — 10.7° found here (Table 1) yields 
the most precise measurements according to our simulations. 


8. Spatial Resolution 

The accuracy of stereoscopic altitude measurements apparently depends on the 
spatial resolution Ax of the instrument, the spacecraft-separation angle a s , and 
the altitude range h max of the solution space, according to the relationship of 
Equation (14), and thus the best spacecraft-separation angle depends on the 
same parameters [Ax] and [/i max [. We determined the minimum of the func¬ 
tion a s (Ax,/i max ) numerically and tabulate the values for Ax = 0.6", 1.6" and 
ftmax = [0.1,0.2]R 0 in Table 1. The coarsest spatial resolution of 2.6" (with pixel 
size of 1.6") corresponds to the EUVI/STEREO instruments, and 1.4" (with 
pixel size of 0.6") to the SDO/AIA instrument with 4 k x 4 k CCD cameras. 

The values in Table 1 show that the spacecraft-separation angle increases 
from Obest = 6.3° to Obest = 10.7° for the altitude range of h max = O.IRq, while 
the accuracy of the stereoscopically triangulated altitudes worsens steadily from 
Ah = 5.5 Mm to 8.9 Mm. Thus, the highest accuracy is clearly achieved for the 
instrument with the highest spatial resolution, as expected. We can derive an 
approximate relationship for the optimum spacecraft-separation angle Obest by 
setting the uncertainty due to the spatial resolution (Ah res ; Equation 11) equal 
to the uncertainty due to loop pairing ambiguities (Ah am b; Equation 12 and 
13). For small spacecraft-separation angles we can then use the approximations 


SOLA: ms.tex; 16 June 2015; 0:33; p. 22 







Blind Solar Stereoscopy 


sin(a s ) ~ a 3 , and cos(ai s ) ~ 1, which yields the relationship 


otbest ~ \ t - > [rad], (15) 

V "■max 

which tells us that the spacecraft-separation angle [a s ] scales with the square 
root of the spatial resolution [Ax]. Varying the altitude range /i max by a factor 
of two, the stereoscopic accuracy worsens a factor of \[2. Thus, for a typical 
range of h max r; (0.1 — O.2)R 0 , and for a pixel size 0.6" ( i.e. Ax = 0.000625) 
as used for the AIA instrument, the optimum spacecraft-separation angle is 
ttbest ~ 4.6° — 6.3°. We repeat the calculations, shown for a spatial resolution of 
1.6" shown in Figure 6, for the AIA spatial resolution of 0.6" in Figure 9. If we 
use the existing STEREO/EUVI instruments with a spatial resolution of 2.6" 
(with pixel size of 1.6", or Ax = 0.0017), the optimum spacecraft-separation 
angle would be in the range of a s ss 7.4° — 10.7° (Table 1). 


9. Spacecraft Position 

Considering the functional behavior of A h(a s ) as shown in Figures 6 and 9, we 
also expect a minimum of the stereoscopic error at large angles of a s = ±90°. 
This solution corresponds to an orthogonal view of an active region, where the 
x and z position of a loop can be measured with maximum accuracy. This large- 
angle configuration, however, has at least three major disadvantages compared 
with the optimum small-angle configuration: (i) Obscuration by the solar limb 
has a much higher probability that an active region is not seen simultaneously 
by two spacecraft; (ii) The stereoscopic correspondence ambiguity is much more 
severe at large spacecraft angles than at small ones; and (iii) A large spacecraft 
distance from Earth demands a much higher telemetry power. For instance, a 
spacecraft-separation angle of a s = 6° implies a ten times smaller proximity to 
Earth than a spacecraft position of a s ~ 60° at Lagrangian points L4 or L5, and 
thus enables a 100 times higher telemetry rate, or a 100 times smaller telemetry 
power for the same data rate. 


10. Non-Potentiality of Magnetic Field 

The ultimate goal of a new stereoscopy mission is a reliable and accurate method 
to measure the magnetic field of solar active regions in order to monitor the 
evolution of non-potential fields and to determine their free energy that can be 
released in large solar flares and coronal mass ejections. The question arises then, 
whether the proposed automated tools are sufficiently accurate for this task. We 
may ask whether a stereoscopy method is sufficiently sensitive to distinguish 
between a potential and non-potential magnetic field, and how accurately can 
it measure the degree of non-potentiality. For a qualitative demonstration we 
simulate the magnetic field for a potential field, using 100 unipolar magnetic 
charges to represent the field (Figure 10; top left), and simulate a corresponding 
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Figure 10. Automated loop tracing with the OCCULT-2 code (red curves), superimposed 
on the simulated magnetic field lines (blue curves in top panels) and the simulated EUV 
images (bottom panels). The left panels show a potential field, while the right panels show a 
non-potential field. Note that the tracings of loops in the EUV images (bottom) are significantly 
different for the potential and non-potential fields, which proves that our automated method 
is very sensitive to the degree of the magnetic non-potentiality. 


EUV image in the same way as described above (Figure 10; bottom left), to 
which then our OCCULT-2 code is applied to trace the loops in an automated 
way (red curves in Figure 10, bottom left). 

Then we use the same model with 100 unipolar magnetic charges, but add a 
vertical electric current to the strongest magnetic charge, which has a magnetic- 
field strength of -2070 G and is located at ( x,y,z ) = (-0.13,-0.26,0.95) (see 
left-most sunspot in Figure 10) by assigning a force-free parameter a-parameter 
of an = 2nN iw i ai /L = 10 -9 cm^ 1 . The method of parameterizing the nonlinear 
force-free field is identical to that of the parameterization used in the simulation 
of the magnetic-field data described in section 3 (Aschwanden et al., 2012a). 
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Comparing this non-potential field (Figure 10; top right) with the potential 
field (Figure 10; top left) one can see that the non-zero a-parameter induces a 
helical twist in the eastern-most sunspot, which causes the field lines emanating 
from the sunspot to be rotated by some amount in the anti-clockwise direction. 
The bottom panels in Figure 10 clearly demonstrate that the automated loop 
tracing code picks out significantly different geometries for the two cases, so that 
we expect a significantly different magnetic field solution, once 3D nonlinear 
force-free modeling is attempted. 

Using the vertical-current model for non-potential magnetic fields (i.e., Priest 
1982; 2014; Aschwanden 2013a), the azimuthal angle \fj\ of a helically twisted 
magnetic field line scales proportionally to the number iVtwist of twists and 
reciprocally to the length L of the magnetic-field line, 

tan(/x) = ^RNtunst _ , (16) 

where R is the flux-tube radius, and a is the nonlinear force-free parameter. From 
the examples of loop tracing using our OCCULT-2 code, as shown in Figure 2, 
we can estimate that the deviation of tracing from a true magnetic-field line 
is less than about one pixel for a loop length of more than 100 pixels. Thus, 
from tan (/ u ) < 0.01 and R w one pixel, we find a lower limit of the force-free 
parameter a = tan (n)/R ~ 1 x 10 -10 cm -1 . This is an estimate of the sensitivity 
of our automated loop-tracing code to the non-potentiality of the magnetic field. 
This is commensurable with a range of a-values found in flare-prone regions, e.g. 
\a\ < 0.02 arcsec -1 ss 3x 10~ 10 cm -1 ; see Figure 3 in Malanushenko et ai, 2014). 

11. Discussion and Conclusions 

We developed an automated stereoscopy code that reconstructs the 3D geometry 
of coronal loops in a solar active region based on EUV images, observed with 
two or three spacecraft located in feasible orbits, i.e. at about 1 AU from the 
Sun, either ahead or trailing behind Earth. In principle, such a “blind stere¬ 
oscopy” algorithm can be applied to already existing spacecraft data, such as 
from SDO/AIA and STEREO/EUVI-A(head) and B(ehind), but these existing 
spacecraft missions have neither an optimum geometric configuration nor opti¬ 
mum spatial resolution of the instruments. We therefore explored the optimum 
conditions for a future mission that may place dual spacecraft anywhere between 
a near-Earth position and the Lagrangian L4 and L5 points. From our study we 
arrive at the following conclusions: 

i) Spatial Resolution : The accuracy of stereoscopic triangulation scales directly 
proportional to the spatial resolution of the EUV imager for a given space- 
craft-separation angle, i.e. the uncertainty or error A h les oc As (Equation 
11). It is therefore desirable to have the highest possible spatial resolution. 
A combination of requirements of full-Sun coverage, high S/N, telemetry, 
and available space-qualified CCD detectors, plus an existing EUV Earth- 
perspective imager already on orbit suggests that a pixel size of 0.6 ,/ for 
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a 4k x 4k imager as implemented in the current SDO/AIA instrument is 
suitable for the purpose of stereoscopic loop tracing. 

ii) Number of Spacecraft: Minimum stereoscopy can be performed with a pair 

of two spacecraft, for instance a spacecraft with one AIA-likc telescope at a 
stereoscopic vantage point, in combination with the existing AIA spacecraft 
in a near-Earth orbit. However, a three-spacecraft configuration, such as the 
existing AIA and two twin spacecraft located ahead of and behind the Earth 
would substantially reduce ambiguities in the stereoscopic correspondence 
problem, as well as provide redundancy in case of any one instrument failing. 

iii) Stereoscopic Correspondence or Ambiguity Problem: From the simulations 
of synthetic images of active regions we found that the number of stereo- 
scopically corresponding loops (detected in an image pair from a spacecraft 
A and B) increases linearly with the spacecraft-separation angle a s . This 
means that larger separation angles lead to increased mapping ambiguity. 
Thus we should aim for the smallest possible angle where stereoscopy is 
feasible. In contrast, large-angle stereoscopy, although feasible, is not at 
optimum conditions. The minimization of the spacecraft-separation angle 
implies also an optimum telemetry rate, since the signal weakens with the 
squared distance to Earth. 

iv) Optimum Spacecraft Separation Angle: The number of ambiguous loops 
and thus the uncertainty of stereoscopic triangulated altitudes increases 
roughly linearly with the separation angle, i.e. A/i am b oc sin(|a a |) (Equation 
12) for small angles. On the other hand, the uncertainty of stereoscopic 
triangulated altitudes decreases reciprocally with the spacecraft-separation 
angle due to the limited spatial resolution of the instrument, i.e. A h les oc 
l/sin\a s \ (Equation 11). The best compromise between these two compet¬ 
ing effects is at a spacecraft separation angle «best where the two uncer¬ 
tainties are comparable, i.e. a.b es t ~ \/2Aa;//i m ax (Equation 15), which 
yields Pbest = 4.6° — 6.3°, for an altitude limit of /i max = 0.1 — 0.2 Rq. 
The altitude limit defines the vertical extent of the 3D solution space in 
stereoscopic triangulations. 

v) Magnetic Non-Potentiality: We performed a qualitative test of how sensitive 

our automated stereoscopy is to the magnetic topology of potential and 
non-potential fields and found that the uncertainties of automated loop 
tracing (in the xy-plane) has sub-pixel accuracy, while the displacements of 
loops between a magnetic potential and non-potential field is much larger 
(of order < 0.1 -Rq), and thus our automated stereoscopy code is sufficiently 
sensitive to changes in the nonlinear force-free field geometry, down to a 
nonlinear force-free parameter of a > l(D 10 cm” 1 . 

Our study complements and augments other recent mission concepts that 
are proposed as platforms to support space weather research and operations. 
One such mission is an L5 Lagrangian Point capability (Vourlidas, 2015) which 
focuses on the propagation of CMEs through the high corona and into the solar 
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wind out to Earth, focusing primarily on the observation and eventual modeling 
of the time and velocity of arrivals of CMEs at Earth. The goals of such an L5 
mission would clearly benefit from a better specification of the Sun-heliosphere 
field interface, for instance by adding magnetograph capabilities to increase the 
coverage of the solar surface. At this point, the proposed concept recognizes that 
’’the entrained magnetic field of an Earth-directed CME is beyond the reach of 
current remote-sensing capabilities.” 

The goal of the present article is to open up a pathway that addresses that 
problem by developing, at least in principle, the methodology to obtain the 
information needed for active-region field modeling based on stereoscopic mea¬ 
surements. Another proposed mission concept, OSCAR (Strugarek et al, 2015), 
for example, also explores that through a perspective from somewhere around 
L5. OSCAR’s premise is that a suitable stereoscopy angle lies in the range of 
« 22° — 125°, based on an earlier work by Aschwanden et al., (2012b). 

The earlier study by Aschwanden et al, (2012b) developed a “quality” metric 
for stereoscopy, which included a plausibility argument for the ambiguity of 
loop tracing. Here, we quantify a metric for ambiguity directly from the loop- 
correspondence algorithm tested here, based on simulated coronal images of an 
active region from different perspectives. A result of the blind-stereoscopy algo¬ 
rithm developed here is that the best performance in resolving the stereoscopic 
correspondence ambiguity is found for considerably smaller spacecraft-separation 
angles of ft 5°. At such small angles, the highest accuracy is found for coronal 
stereoscopy to constrain nonlinear force-free modeling of coronal magnetic fields. 
This low separation angle, additionally offers the advantage of a higher telemetry 
rate for a given spacecraft and ground-antenna combination. 

With these findings we realize that any space-weather mission concept has 
to deal with the trade-offs between small-angle stereoscopy, large-angle CME 
coverage, and spacecraft telemetry. An optimum configuration that meets the 
scientific and operational needs to study CMEs both in terms of magnetic con¬ 
tent and in terms of their propagation through the heliosphere (c/., Schrijver et 
al, 2015) would appear to substantially benefit from having a triplet spacecraft 
configuration that includes both the LI and L5 point, and a third spacecraft 
about 5° away from Earth, equipped with complementary capabilities to measure 
the magnetic field, the plasma parameters, the arrival velocity, and the arrival 
time of Earth-bound (geoeffective) CMEs. 
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